Vortex induced deformation of the superconductor crystal lattice 
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Deformation of the superconductor crystal lattice caused by Abrikosov vortices is formulated as a 
response of the elastic crystal lattice to electrostatic forces. It is shown that the lattice compression 
is linearly proportional to the electrostatic potential known as the Bernoulli potential. Eventual 
consequences of the crystal lattice deformation on the effective vortex mass are discussed. 
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During the transition from a normal to a supercon- 
ducting state, metals reduce their specific volumes 0,0- 
In a mixed state, would it be the Abrikosov vortex lat- 
tice or a structure of lamellas, the superconductivity is 
locally suppressed and the specific volume is inhomoge- 
neous. The mixed state is thus accompanied by strains 
and stresses, which enter the balance of total energy. 

In general, the energy of strains is much smaller than 
the energy of the superconducting condensation and the 
energy of the magnetic field. Its contribution becomes 
appreciable only under special conditions. For example, 
experiments on single crystals of Pb-alloys 0, 0] and Nb- 
alloys [J, la, la] revealed that an orientation of the vortex 
lattice is influenced by its angles to main crystal axes. 
Since the gap of alloyed samples is quite isotropic, purely 
electronic models have failed and this effect has been ex- 
plained with the help of strains induced by vortices Q. 

In the 1990th a different structural effect has been ob- 
served on NbSc 2 . If the magnetic field is tilted from 
the c-axes, the Ginzburg-Landau (GL) theory predicts a 
state in which rows of vortices are aligned with the direc- 
tion of tilting , while experiments 0, 0, ^| show 
them aligned in the perpendicular direction. Again, the 
interaction of vortices with the crystal strain explains the 
observed alignment [liij . 

Finally, we would like to mention phenomena which 
are predicted but not yet fully experimentally confirmed. 
Perhaps, one of the most interesting predictions is a siz- 
able contribution of the lattice deformation to the mass 
of vortex [3 13 EH 13- Besides, there is a number of 
phenomena due to strains at surfaces which are discussed 
in Ref. E3- It is a l so argued that the strain can medi- 
ate an attractive long-range interaction between vortices 

0- 

As far as we know, all theoretical studies of deformable 
superconductors use a phenomenological model, which 
assumes that the superconducting condensate interacts 
with the lattice density directly via a strain dependence 
of material parameters . This model dates back to the 
1960th [H H3> when it 

was used to describe the vortex 
pinning. The strength of the interaction is deduced from 



changes of the specific volume in the phase transition, see 
also Ref. Q. 

In this paper we assume that the condensate interacts 
with the crystal lattice via electrostatic forces created 
by the so called Bernoulli potential. We show that this 
mechanism results in the interaction based on the specific 
volume. In addition to known theories we obtain gradi- 
ent corrections and demonstrate that they are important 
for the vortex motion of the Abrikosov vortex lattice in 
niobium. 

In an isotropic continuum, the displacement field u 
obeys the equation 22] 
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H V(V.u) -(iVxVxu = F, 



(1) 



where K and \i are the bulk and shear modulus, and F 
is the volume density of force acting on the lattice. Some 
authors prefer to express coefficients on the left hand 
side in terms of the Poisson ratio a — (3K — 2fx)/(6K + 
2fi) and the Yang modulus E = 3K(l — 2a). In the 
basic approximation K and \i are constants. Their small 
change in the superconducting state was assumed in |23| . 

The inhomogeneous superconductivity results in the 
force F. The present theory differs from previous ap- 
proaches in the approximation adopted for F. Let us 
sketch the approach based on the specific volume first. 
More details the reader can find in Ref. Ell- 
in the phase transition from the normal to the super- 
conducting state, a system shrinks by a volume difference 
SV = V n — V s = a? V. A typical value of or is about 
10~ 7 . The density of the atomic lattice correspondingly 
increases, 5ni a t = ^f at — ^" at = cut n. 

The strain coefficient ut depends on the temperature 
via a fraction w of electrons, which become supercon- 
ducting, ax = ato. Here a is the strain coefficient at 
zero temperature. In the spirit of the GL theory, we ex- 
press the superconducting fraction in terms of the GL 
function id = 2\ip\ 2 /n, i.e., Sni at — 2a \ip\ 2 . 

In the vortex core or in the region of surface currents, 
the GL wave function changes in space, see Fig. ^ The 
lattice then tends to be inhomogeneous, which causes 
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internal stresses. These stresses lead to a density of force 

F Sim = KaV^^ (2) 
n 

proposed by Duan and Simanek [lfij in their study of 
the vortex mass. The phenomenological formula fur- 
nishes us with a density of forces with no regards to how 
the superconducting electrons are coupled to the lattice. 

Let us try to express such force on the lattice in a semi- 
microscopic way. Diamagnetic currents, either being on 
the surface or circulating around the vortex core, always 
cause inertial and Lorentz forces, which are balanced by 
an electrostatic field E = — V</>, see London 24]. This 
electric field transfers the Lorentz force from electrons to 
the lattice, therefore one can expect that it also causes 
lattice deformations. Accordingly, we suppose that the 
electrostatic field force 

F = enV0, (3) 

is playing the role of the force F in equation 0J. For 
simplicity of notation we assume singly ionized atoms, 
e is the charge of an electron, so that the ionic charge 
density is — en. 

The electrostatic potential <j> is known as the Bernoulli 
potential. It has been derived in a number of approxi- 
mations [2 4 . l25l |26L |27| . Here we will use the formula of 
Ref. [13 

1 2 
efi = ~ o * ip(-ihV~e*A) tp 
2m*n 

The space profile of the potential cj> is shown in Fig. [3] 
and the individual terms to the potential are compared in 
Fig. El The first term in (@J is the quantum kinetic energy 
and represents the gradient corrections. In the London 
limit it reaches the form of the classical Bernoulli law, 
4>(-ihV - e*A) 2 i>/2m*n -> e* 2 A 2 to/4m* = Lumv 2 /2, 
which gave the name to the entire potential. The su- 
perconducting fraction lo multiplying the kinetic energy 
accounts for the fact that the Lorentz and inertial forces 
act exclusively on the super-electrons, while the balanc- 
ing electrostatic force acts on all electrons 25] . This 
force being proportional to the square of the magnetic 
field has been used to calculate the shape distortion by 
flux-pinning-induced magnetostriction [28j . 

The second term of the Bernoulli potential (@J is the 
dominant one and we will focus our discussion on it. 
This second term is identical to the potential derived 
by Khomskii and Kusmartsev [2(| from the effect of the 
BCS gap on the local density of electronic states. The 
third term of 0} we call the entropic correction. 

It is worth noting that all the three equations ® 
and Q are resulting from Gibbs variational principle if 




FIG. 1: The superconducting fraction uj — 2\ip\' 2 /n in the 
triangular Abrikosov vortex lattice. We assume niobium with 
the GL parameter increased by non-magnetic impurities to 
ft = 1.5, the temperature T = 0.7 T c and the mean magnetic 
field B = 0.24 B C 2- In centers of vortices the superconducting 
fraction to goes to zero, while between vortices it approaches 
its non-magnetic value 1 — T 4 /T c 4 = 0.76. 



the free energy used in [2jj is extended by the ionic lattice 
deformation energy. 

Let us link the deformation caused by the electrostatic 
field with the standard theory of magnetostriction. At 
zero temperature, the Gibbs energy of normal and su- 
perconducting states differ by the condensation energy 
0, G s = G n - Ve con , where e con = 7 T c 2 /4 = B 2 /2n . 
Since the pressure derivative of the Gibbs energy deter- 
mines the sample volume, V sn = dG Stn /dp, one finds 
V s = V n — V de con /dp, i.e., a = de con /dp. This relation 
allows us to show that the force Fsi m from © equals 
the electrostatic force caused by the second term of the 




FIG. 2: The electrostatic potential for the parameters of 
Fig-0 Its shape reminds the superconducting fraction shown 
in Fig.0 which indicates that corrections beyond the approx- 
imation of Khomskii and Kusmartsev are small. 
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FIG. 3: Cuts through the electrostatic potential from Fig. [5] 
The thick full line represents the total potential JIJ, the dot 
line is the Khomskii-Kusmartsev approximation. The thin full 
line is the gradient correction given by the 'kinetic energy' 
term of and the dashed line is the entropic correction 
given by the last term of (J1J. 



Bernoulli potential Q. 

To proceed, we will use the fact that the pressure mod- 
ifies the condensation energy indirectly by an increase of 
the electron density 



de r 



de rm dn de c . 



dp 



dn dp 



dn K 



(5) 



In the rearrangement we have employed the definition 
of the bulk modulus K ~ —V/(dV/dp) = n/(dn/dp). 
Substituting (0 into the force J5J) one finds 



Sim 



de CO n y 2h/f 
dn n 



(6) 



Comparing © with J3J one can see that the force in- 
troduced by Simanek equals to the electrostatic force due 
to the second term of the Bernoulli potential 0}. In 
this sense, the approximation used within the theory of 
dcformable superconductors is equivalent to the approxi- 
mation of the electrostatic potential derived by Khomskii 
and Kusmartsev. The gradient and entropic corrections 
of the Bernoulli potential 0} provide us with correspond- 
ing corrections to the force of Simanek. 

More interesting is the gradient correction. From the 
formula follows that far from the vortex core, where 
ipip — > (1 — T 4 /T c 4 ) n/2, the gradient correction is propor- 
tional to the square of the local current. For an isolated 
vortex, the gradient correction thus decays on the scale 
of the London penetration depth. Consequently one can 
expect, that it plays an important role in high k materi- 
als. 

The effect of the gradient correction is traceable also 
for a conventional material assumed here. To be specific, 
our sample is a niobium rod parallel to the magnetic 



field. The GL coherence length of niobium is reduced 
by non-magnetic impurities so that the GL parameter is 
increased to k — 1.5, while other material parameters re- 
main close to values of the pure niobium. All plots are 
for T = 0.7 T c and B = 0.24 B c2 . In this case the mag- 
netic field is not split into separated unitary fluxes but 
it is nearly homogeneous with amplitude fluctuations of 
about 20% around the mean field B. The vortex cores are 
well separated, however, as lu reaches its non-magnetic 
value in the out-of-core region, see Fig. ^ 

Let us analyze the three contributions to the electro- 
static potential from the point of view of forces they 
cause. According to the position of the inflex point of 
the Khomskii-Kusmartsev potential seen in Fig. |3 one 
can estimate that the maximum of the Simanek force is 
at about x,y ~ 0.1a. This is quite close to the cen- 
ter of the vortex core. The gradient correction oscillates 
fast in space having a magnitude much smaller than the 
Khomskii-Kusmartsev potential. Its gradient, however, 
is rather comparable to the gradient of the dominant 
term, which shows that gradient contribution to the force 
can appreciably modify the Simanek force. In the heart 
of the vortex core, the gradient correction to the force 
acts against the Simanek force, while in the skin of the 
vortex core it points in the same direction. As a result, 
the maximum of the total force is shifted outwards to 
x,y ~ 0.2 a. 

Naturally, gradient corrections modify the vortex 
mass. Figure 0^ shows how individual regions in the 
Abrikosov vortex lattice contribute to the vortex mass 
for the Simanek force, in figure 0Jd the gradient and the 
entropic corrections are included. The plotted function 
is the density of kinetic energy of lattice ions driven by 
vortices moving with velocity V in the x direction [bH ] 



E kin = ^V 2 nM 



du x 
dx 



du y 
dx 



(J) 



where M is the mass of a single ion. 

In most places of the Abrikosov vortex lattice, the ionic 
kinetic energy is lowered by the correction terms, see 
Fig. 01 From the integral over the elementary cell one 
obtains the vortex mass per unit length. Keeping both 
corrections, the ion contribution to the vortex mass is 
reduced by a factor 0.83 as compared to Simanek result. 
If one neglects the gradient correction keeping only the 
entropic correction, the reduction factor is 0.70. The gra- 
dient correction thus leads to a small enhancement of the 
vortex mass. 

Comparing relative amplitudes of the kinetic energy 
inside and between the cores, one can see that the correc- 
tions have increased the share of the out-of-core region. 
Since the entropy term merely reduces the amplitude of 
the Simanek force, this redistribution of distortions is 
exclusively due to the gradient correction. According to 
Cano, Levanyuk and Minyukov 19] the out-of-core lat- 
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FIG. 4: The density of kinetic energy of lattice ions created 
by vortices moving in the x direction for parameters of Fig. Q 
The kinetic energy due to (a) Simanek force © , and (b) the 
electrostatic force The plotted function is dimensionless 
a quantity identical to the bracket from the formula (7|. 



tice deformations are important for the strain-mediated 
interaction of vortices. We suggest that the theory of 
this interaction should be reexamined with the gradient 
correction included. 

In summary, we have expressed the forces deforming a 
lattice of a superconductor in terms of the electrostatic 
force. This approach allows one to benefit from expe- 
riences accumulated within the theory of the so called 
Bernoulli potential. One directly obtains the gradient 
corrections we have discussed in this paper. We suggest 
that the theory of vortex motion should be revised tak- 
ing this gradient corrections into account. The present 
theory is restricted to homogeneous isotropic materials. 
Its extension to layered materials we plan to discuss in 
short future. 
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